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Ytterbium divalency and lattice disorder in near-zero thermal expansion YbGaGe 
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While near-zero thermal expansion (NZTE) in YbGaGe is sensitive to stoichiometry and defect 
concentration, the NZTE mechanism remains elusive. We present x-ray absorption spectra that 
show unequivocally that Yb is nearly divalent in YbGaGe and the valence does not change with 
temperature or with nominally 1% B or 5% C impurities, ruling out a valence-fluctuation mechanism. 
Moreover, substantial changes occur in the local structure around Yb with B and G inclusion. 
Together with inelastic neutron scattering measurements, these data indicate a strong tendency for 
the lattice to disorder, providing a possible explanation for NZTE in YbGaGe. 

PACS numbers; 65.40.De, 61.10.Ht, 75.20.Hr, 63.20.-e 
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Observations of near-zero thermal expansion (NZTE) 
are relatively rare, especially near room temperature. 
The most famous example is invar (Fe64Ni36) which has 
a volume expansion cocfhcient (3 « 4 x 10^® K^^. Re- 
cently, YbGaGe (Fig. [1]) has been identified as a poten- 
tial NZTE material,! with /? w -3 x 10"'^ K'^^ depend- 
ing on the Ga/Ge ratio. In fact, the exact stoichiometry 
is an important factor, with many subsequent measure- 
ments unable to reproduce the original worki^i^i^i^ii In 
particular, (3 varies widely between about 4.5 x 10~^ K^^ 
and 1.4 X 10"'* K"*, typical values for intermetallics. Un- 
derstanding the differences between these samples should 
eventually allow the reproducible fabrication of, what 
should be, a technologically important material. 

A fundamental difference between the original and sub- 
sequent measurements is the magnitude and the sign of 
the magnetic susceptibility. In Ref. J, a Curie- Weiss-like 
susceptibility was observed above about 100 K with an 
effective moment of fi^s ~ 4.12/iB, close to the full mo- 
ment expected from free Yb'^+ ions (4.54 /^b)- Below 
100 K, a sharp decline in the moment occurred, with 
fi^ff « 0.82/zb. These authors pointed out that a bond- 
valence sum® for each of the two ytterbium sites (Fig. [1]) 
gives a valence of -1-2.6 for the Yb(l) site and -1-2.0 for the 
Yb(2) site. Such a result is consistent with a mixed va- 
lent state for at least one of the Yb sites. Together with 
the magnetic susceptibility data, this observation led the 
authors to conclude that the mechanism for the observed 
NZTE was a change in the Yb valence with temperature 
toward a divalent state below 100 K, consistent with the 
susceptibility. 

In direct contrast to those first measurements, subse- 
quent measurements find a diamagnetic susceptibility for 
the pure compound (^i^i^i^ even over a wide range of the 
Ga/Ge occupancy ratic^ and with other defects^ These 
susceptibility measurements are a very strong indication 
of divalent Yb at all measured temperatures. However, if 
the real Yb valence is -1-2.6 or even less, this could be an 
indication of a very high Kondo temperature, Tk, per- 



haps well in excess of 1000 K. Given this possibility and 
the fact that the magnetic susceptibility is expected to go 
to a constant xo oc l/Tk at temperatures below Tk, the 
diamagnetic contribution from the nonmagnetic matrix 
needs to be carefully considered before declaring that all 
the Yb is divalent. A spectroscopic measurement of the 
Yb valence would provide a direct measurement. 

As a possible explanation for the apparently anoma- 
lous Salvador et al. work, Drymiotis et al7 considered 
the possibility that the observed NZTE could be a conse- 
quence of disorder. They note that the original synthesis 
was performed in graphite crucibles and therefore grew 
YbGaGe with small amounts of carbon and boron in alu- 
mina crucibles. They found that the thermal expansion 
could be reduced by about 50% with only 0.5 at.% car- 
bon or boron included in the starting materials. (All 
subsequent concentration percentages are in at.%). They 
conjecture that, although they did not obtain NZTE and 
still observe diamagnetic behavior, that the mechanism 
for NZTE in Ref. 1 has a similar origin. 

Here, we report x-ray absorption near-edge structure 
(XANES) measurements at the Yb Lni edge that show 
ytterbium in these materials is, in fact, nearly divalent, 
that the valence does not change dramatically with car- 
bon or boron doping, and that the remaining trivalent 
component is easily explained as due to an impurity 
phase (for example, Yb203 or a Yb-Ga binary alloy). In 
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FIG. 1: Stacking of the Yb(l)-Ga6 trigonal prisms and the 
Yb(2)-Ge6 octahedra in the PQ^/mmc hexagonal cell. 
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FIG. 2; Yb Liii-edge XANES spectra for the pure, 1% B, 
and the 5% C samples of YbGaGe, all collected at 30 K, 
together with a room-temperature spectrum of Yb2 03. The 
three YbGaGe spectra are nearly identical, and are therefore 
difficult to discern in this figure. 



addition, extended x-ray absorption fine-structure (EX- 
AFS) measurements of the local structure environment 
around the Yb atoms shows that small amounts of car- 
bon or boron in YbGaGe have a surprisingly large effect 
on the average local structure. This extreme sensitivity 
to interstitials appears to be the root cause of near-zero 
thermal expansion in this system. 

Samples were prepared as in Ref. 7. For the present 
work, in addition to a high purity sample grown in a yt- 
tria crucible, nominally 1% boron and 5% carbon doped 
samples were also fabricated. X-ray diffraction measure- 
ments indicate the YPtAs structure type in the hexago- 
nal PQs/mmc space group, with small (<5%) uniden- 
tified impurity phases. The yttria grown sample has 
P « 4.6(1) X 10~^ at room temperature accord- 

ing to powder diffraction measurements.^ We have re- 
arcmelted the doped samples to obtain a polycrystalline 
average and obtained /3 = 3.3(3) x 10"^ K'^ for the 1% B 
sample and 4.4(2) x 10"^ for the 5% carbon sample 
using a capacitance dilatometer.— These values are com- 
parable to the values obtained with diffraction for the 
0.5% B and C doped samples (3.0(3) x IQ-^ and 
3.2(1) X 10"^ K-\ respectively^). Although the doped 
samples have been re-arcmelted and are compared us- 
ing different measuring techniques, these data counter 
the notion that more boron or carbon would decrease (3. 
Magnetic susceptibility measurements indicate diamag- 
netic behavior for all samples with a small impurity tail 
corresponding to free Yb^^ moments. The effective mo- 
ment per mole and an estimate of the crystal field effects 
indicates the Yb^+ impurity fraction is between 1-5% and 
shows no correlation with B or C concentrationi^i 

XANES data were collected on Beamline 10-2 at 
the Stanford Synchrotron Radiation Laboratory (SSRL). 
The samples were first re-ground, passed through a 20 



/im sieve, brushed onto adhesive tape, and stacked to 
achieve an absorption length change at the Yb Lni edge 
(8.944 eV) of A/ii « 1. The samples were then placed 
in a LHc flow cryostat and measured between 30 and 
300 K. Data were collected in transmission mode with 
a defocused beam from a half-tuned Si(220) double- 
crystal monochromator resulting in an energy resolution 
of about 2.0 eV. Since the core-hole lifetime is about 4.2 
eV, the reported spectra are not lifetime limited. All 
spectra were collected under the same conditions. Data 
were reduced and fit using standard procedures. In par- 
ticular, the absorption from the Yb Lni edge as a func- 
tion of the incident energy E, iJig_{E), was isolated from 
the total absorption ijl{E) after approximating all other 
contributions by extrapolating from the pre-edge data 
and forcing /ia(i^) above the edge to follow a Victoreen 
formula.— The EXAFS oscillations were isolated by fit- 
ting a 7-knot cubic spline to approximate the embedded- 
atom absorption /ig {E) , resulting in the EXAFS function 
x{k) — ifJ-aik) — ^o{k))/^o{k), where k is the photoelec- 
tron wave vector obtained from E and the threshold en- 
ergy Eq, arbitrarily taken as the energy at the half- height 
of the edge. The oscillations in x{k) are due to interfer- 
ence of the back scattered and outgoing parts of the pho- 
toelectron wave function modulating the absorption coef- 
ficient. A Fourier transform (FT) of x{k) produces peaks 
in r-space corresponding to scattering shells around the 
absorbing atomic species. Note that phase shifts of the 
outgoing and back scattered photoelectron make the FT 
more complicated than a true radial distribution func- 
tion, so detailed fits using the RSXAP packageiSiiiii^ 
and theoretical lineshapes calculated by FEFF7— were 
performed to extract the local structure information. 

The XANES results are displayed in Fig. [2l together 
with the spectra of Yb203 as a trivalent Yb reference. 
The main peak ("white line") position of the Yb203 ref- 
erence is typical of trivalent Yb, including in intermetallic 
compounds. There is a clear energy shift of about 8 eV 
between the YbGaGe compounds and Yb2 03, indicating 
Yb in these intermetallics is predominantly divalent. In 
addition, we observe a small bump at the Yb2 03 white 
line position, indicating a small trivalent component. By 
fitting these features to a pseudo-Voigt function each for 
the divalent and the trivalent resonances and an inte- 
grated pseudo-Voigt to model the edge step, we obtain 
the /-hole occupancy rif = 0.17 ±0.10. The small triva- 
lent component could be due to an intermediate valent 
ground state in YbGaGe, but is more likely due to in- 
termetallic or oxide impurities. Very little difference in 
the /-occupancy is observed between the three samples. 
Moreover, there is very little change in with temper- 
ature, less than 1% between 30 and 300 K. 

The EXAFS data and fit for the pure YbGaGe sample 
are shown in Fig. [3)^ The data quality is very high and 
the r-space fit quality is also very good above about 2 A. 
However, below 2 A, there is some unexplained ampli- 
tude, as reflected in the reported R{%) = 13.1 (a high- 
quality EXAFS fit should be near 5%). Some of this 
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FIG. 3: EXAFS data in r-space for (a) the pure YbGaGe 
sample, and (b) all three samples at 30 K. The Fourier trans- 
forms (FT) are between 2.5 and 16.0 A~^, Gaussian narrowed 
by 0.3 A"^ The fit in (a) is between 1.2 and 6.0 Aand the 
outer envelope is ± the amplitude, and the modulating inner 
line gives the real part of the complex transform. Only the 
amplitude is shown in (b). 



amplitude can be explained with between 5-10% Yb203. 
Even so, the fit shown in Fig. [3lja) aheady includes such 
an impurity, so there is likely some other kind of impurity 
present, possibly in addition to the sesquioxide. We have 
tried including other intermetallic phases in the fit, such 
as YbGa and YbGa2, with only moderate success. In any 
case, the fit parameters are very close to the expected val- 
ues from bulk diffraction studies, and the variance of the 
pair-distance distribution widths (Debye- Waller factors) , 
cr^, are also all of a reasonable magnitude. 

The r-space transforms for all samples at 30 K are 
shown in Fig. ^h). There are clear differences between 
the pure and the substituted samples at all length scales. 
The fit results^ for the substituted samples are similar to 
the pure case with the following exceptions: (1) the bond 
lengths are very similar except for the Yb(l)-Ge pair in 
the 1% boron substituted sample, nominally at 3.35 A, 
but instead at 3.13(3) A; (2) many of the Debye-Waller 
factors are substantially different; and (3) the amplitude 
reduction factor Sq is unphysically small, « 0.5-0.6 in 
both cases. These results point strongly to the pres- 
ence of multiple phases over and above that observed for 
the pure YbGaGe sample. For instance, the short value 
for the Yb(l)-Ge pair may be indicative of a YbGa-like 
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FIG. 4: Neutron scattering intensity versus energy transfer 
AE at 15 K and 300 K for undoped and C-doped YbGaGe. 
The 78 K data are not shown for clarity. The incident energy 
in the main panel was 25 meV, and 125 meV in the inset. 

phaseM- More directly, the reduction in Sq (an overall 
scale factor) and the variations in the Debye-Waller fac- 
tors, while still allowing for a high quality fit, can be 
explained by allowing for nearly amorphous or highly dis- 
torted regions of the sample. The presence of such a high 
fraction of sample being so disordered (about 20%) indi- 
cates that the presence of boron or carbon in the lattice 
has relatively long-range effects on the local structure. 
This result could be verified through quantitative analy- 
sis of x-ray or neutron diffraction data. 

In order to determine whether these changes are vi- 
brational in origin, inelastic neutron scattering measure- 
ments were performed on polycrystalline samples of Yb- 
GaGe, both undoped and doped with 1% carbon. Data 
were collected on the PHAROS time-of-flight spectrom- 
eter at the Manuel Lujan Los Alamos Neutron Sci- 
ence Center (LANSCE) with incident energies of 25 and 
125 meV and at T =15, 78 and 300 K. Groups of detec- 
tors were summed to improve statistics. Figure S] shows 
data from detectors with an average scattering angle of 
= 89°; for such large scattering angles, the spectrum of a 
polycrystalline sample approximates the phonon density- 
of-states. Peaks are observed near 6, 8.5, 16 and 30 meV. 
The fact that the phonon spectra for both pure and 
carbon-doped YbGaGe samples are identical means that 
the changes in the thermal expansion and EXAFS data 
from the C-doped samples cannot arise from a bulk dis- 
tortion of the lattice, since this should affect the phonon 
frequencies. Instead, these changes must be due to ran- 
dom or near-random disorder. Presumably, this disorder 
reduces the anharmonicity in the pair potentials that gen- 
erate the thermal expansion. Apparently, these neutron 
measurements are not sensitive to such an effect. More- 
over, fits to the EXAFS including an anharmonic term, 
that is, a third cumulant,^^ are inconclusive. 

These results strongly support the conclusion by 
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Drymiotis et alX that the mechanism for the reduced 
thermal expansion in B and C doped samples, and possi- 
bly for the NZTE observed by Salvador et al.^ is chemical 
defects. It also spectroscopically confirms the conclusions 
drawn by several researchers^i^ii that the diamagnetic 
susceptibility is indicative of a divalent Yb state at all 
measured temperatures. 

Unexplained facts remain regarding the original Sal- 
vador et alX work. First, no subsequent published mea- 
surement has produced a magnetic susceptibility consis- 
tent with mostly Yb^+ above 100 K. This discrepancy 
may be due to the exact impurities that may have been 
present in the original samples. However, given the mag- 
nitude of the susceptibility, it seems unlikely that an im- 
purity could generate a susceptibility corresponding to 
bulk Yb^+ at room temperature. Given this structure's 
propensity to disorder with small defect concentrations, 
the possibility remains that the original samples did, in 
fact, undergo a valence change upon cooling. 

Another peculiar fact is that the bond-valence sums 
give a much higher valence to the Yb(l) site (-1-2.6) than 
the Yb(2) site (-1-2.0). This analysis depends on the z- 
parameter for the Ga and Ge sites, a quantity that has 
not been re-measured. The EXAFS results reported here 
are similarly insensitive to the differences between the 
Ga and Ge atoms as are the x-ray diffraction measure- 
ments; however, the results are consistent with two sep- 
arate Yb sites with different coordination environments, 
as expected from the diffraction measurements. Indeed, 
the local measured bond lengths are all consistent with 
those expected from the diffraction results and their asso- 
ciated z-parameters. The bond-valence sum conclusion 



also depends on the relative mixture of Ge and Ga on 
their nominal sites, a quantity that is essentially assumed 
both by the diffraction and the EXAFS measurement. 

In conclusion, the /-hole occupancy Uf has been mea- 
sured in YbGaGc using Yb ini-edge XANES spec- 
troscopy, and apart from a small trivalcnt component 
that is likely due to oxide or intermetallic impurities, 
Yb in this system appears to be divalent, consistent 
with several previous measurements, such as magnetic 
susceptibility.— i^i^ii The valence does not change dramat- 
ically either with temperature or with small amounts of 
boron or carbon substitution. Although the local struc- 
ture obtained by fitting the EXAFS is consistent with 
that expected from the bulk crystal structure, systematic 
differences with the fitting model indicate the presence of 
some other phase, especially in the substituted materials. 
In particular, the 1% boron and 5% carbon substituted 
samples show large differences in the EXAFS compared 
to the pure sample, indicating a very large effect on the 
local lattice for such small amounts of impurities. The 
inferred presence of such highly disordered regions is con- 
sistent with the conclusion reached in Ref. that the 
reduction of the thermal expansion is related to defects 
caused by chemical substitution. 
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